Volovik effect on NMR measurements of unconventional superconductors 
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We studied the Volovik effect on the NMR measurements of the two unconventional superconducting (SC) 
states, the d-wave and ±s-wave states. We showed the generic field dependencies of the spin-lattice relaxation 
rate 1/Ti and Knight shift K at low temperature limit in the pure cases as: 1/Ti oc HlogH, K oc VTT for 
the d-wave, and 1/Ti oc H, K oc H for the ±s-wave state, respectively. Performing numerical calculations 
we showed that these generic power laws survive for the good part of low field region with the realistic amount 
of impurities. We also found that the Volovik effect acts as an equivalent pair breaker as the unitary impurity 
scattering, hence induces the same temperature evolutions on 1/Ti and K, respectively, as the unitary impurities 
in both SC states. This finding implies that the Volovik effect should always be taken into account for the analysis 
of the NMR measurements in the mixed state. 

PACS numbers: 74.20.Rp, 74.25.fc,74.25.Uv 



Introduction. - There are various experimental probes of 
the gap symmetry of the superconductors. Mostly, they are 
probing low lying excitations in the superconducting (SC) 
state to distinguish a full gap, a nodal, or a point gap. The 
measurement results usually produce the typical temperature 
dependencies, for example, an exponentially flat for a full gap 
or various power laws for a nodal or point gap. The Volovik 
effect (VE) [1] has enlarged the probing dimension to the 
magnetic field axis. The principle is actually very simple in 
that the low energy density of states (DOS) of the SC state is 
changed in the mixed state with vortices induced by the ap- 
plied magnetic fields and then its field dependence is probed. 
In particular, this change of DOS due to the vortices should 
be more sensitive with a gapless superconductor such as the 
d-wave state in the cuprate superconductors |T| 0] • However, 
recently, it was shown that the strong field dependence of the 
Volovik effect is not unique to the nodal gap superconductors 
but also can appear with the full gap s-wave superconductors 
if there exist multiple s-wave gaps with different gap sizes J^. 

On the other hand, most of the Volovik effect measurements 
have been carried out with the specific heat (SH) and thermal 
conductivity. However, in principle, it should be effective with 
any experimental measurements which probe the low energy 
DOS and its variation with the applied field. Therefore, the 
NMR measurements such as the spin-lattice relaxation rate 
l/Ti(T, H) and Knight shift K(T, H) should also be good 
probes for the Volovik effect. In comparison with the SH 
and thermal conductivity measurements, the NMR measure- 
ments of Volovik effect require much harsher conditions such 
as the sample size, low temperature control, and high qual- 
ity homogeneity of the field, etc. Nevertheless there already 
exist a few experimental works, for example, on Fe pnic- 
tide superconductors, Ba .6gKo.3iFe 2 As 2 by W. R Halperin 
and coworkers J4(] and BaFe 2 As 2 (As .67Po.33)2 by Nakai et. 
a/.j5[], as well as on the high-T c cuprates by G. Q. Zheng and 
coworkers jg]. However, because the theoretical study on this 
subject is very rare |2U, the reliable interpretation and the ex- 
traction of the useful information from the experimental data 
are limited. Therefore more systematic and detailed theoret- 



ical study of the Volovik effects on the NMR measurements 
is highly demanded and it will enhance the capability of the 
NMR technique to study the superconducting gap symmetry. 

In this paper, we specifically studied the NMR 1/T\ and 
Knight shift K in the SC state of the two typical uncon- 
ventional superconductors, namely, the d-wave and ±s-wave 
state slU, Ell- We derived the exact power laws of the generic 
field dependencies of both SC cases in pure state: 1/Ti oc 
H log H,K oc VH for the d-wave, and 1/Ti oc H, K oc H 
for the ±s-wave state, respectively. We also numerically 
studied the impurity effects on 1/Ti and K, using the self- 
consistent T-matrix approximation (SCTA) fioill ill and pro- 
vided the systematic comparison between the d-wave and ±s- 
wave gap states. We found that the impurity scattering sub- 
stantially reduces the low field region where the generic field 
dependencies survive in the ±s-wave case, while it doesn't 
much affect the generic power laws except a constant shifting 
in the d-wave case [12]. Another key finding of our numeri- 
cal study is that the Volovik effect induced by magnetic fields 
practically acts as an equivalent pair breaker as the strong cou- 
pling (unitary limit) impurities, hence modifies the low tem- 
perature behaviors of the 1 jT\ and K in the same fashion as 
due to the formation of the resonant impurity band. This result 
implies that it is always necessary to take the Volovik effect 
into account when one interprets the low temperature NMR 
experiments with magnetic field in order to extract the correct 
information from the data. 

Formalism. - For most of purposes to study the Volovik 
effect, we just need to calculate the position dependent DOS 
N(cj,r) in the presence of vortices. Using the semiclassical 
approximation, the matrix form of the single-particle Green's 
function in the SC state, including Doppler shift of the quasi- 
particle excitations e(k) due to the circulating supercurrent 
v s (r), is given by HQ] 



G(k,r,w) = 



[u + v s (r) • k]r + e(fc)r 3 + An 



v s (r) • k] 2 - e 2 (k) - A 2 



(1) 



where Ti are Pauli matrices and r is the distance from the vor- 
tex core. A is the SC gap function and v s (r) is ~ mr@- 



2 



The position dependent DOS is calculated as N(ui, r) — 
— iTrlm J^k Go (k, r, tS). Finally, the field dependent quan- 
tities are obtained from the areal average DOS per unit vol- 
ume as N(u>,H) = J^ H dr 2 N(uj,r)/-n:R 2 H with the mag- 
netic length Rfi = \J~§fj ($0 a f ux quanta) and the SC co- 
herence length £. 

In the homogenous SC state, the general structure of the 
l/Ti is written as 



and 



and the Knight shift K as 



K = -T 



where N(u) = (Re- 



df FD (uj) 
acu N(u)) 

OLO 



(3) 



is the Fermi surface (FS) 
— \ " is the similar 



V" 2 -A 2 (9) /0 

averaged DOS and M(w) = (Re- . 

quantity induced in the SC state. In the mixed state, we only 
need to replace N(oj) by N(w,H) and M(w) by M(u,H) 
in the above formulas (2) and (3). M(w, H) is obtained from 
M(u! 7 r) = — ^Trlm^ fe Gi(k, r, u) and its areal average. 
Then we can calculate the field as well as temperature depen- 
dent NMR properties l/Ti(T, H) and K(T, H). 

In the mixed state of the d-wave superconductor, it is well 
known that the N(lu — 0,r) ~ 1/r due to the linear DOS 
N(ui) ~ to and the Doppler shifting energy v s (r) • kjr ~ 1/ r. 
Also knowing that M(lo) vanishes due to the gap symmetry, 
we can readily extract the field dependence of the l/Ti for the 
d-wave state from Eq.(2) as 



-L(T^0,/T)~ f " dr 2 (-) 2 /nR 2 H ^ HlogH, (4) 
Ti J$ r 

and similarly for Knight shift as 



dr 2 {-)/nR 2 H 



H. 



(5) 



In the case of the ±s-wave state with small gap A5 and 
large gap A^, again all we need to know is the DOS near 
zero energy in the mixed state. The author[3] has recently 
shown that the Doppler shifting energy v s (r) • kp ~ A^/r 
overshoots the small gap to induce a constant DOS - which is 
basically the normal state DOS of the small gap band Ng orm - 
around ui = in the small gap band in the region around the 
vortex core to the distance r* — f^- that is a field inde- 

"» As 

pendent constant. Hence the field dependencies of the NMR 
properties of the ±s-wave state are the following. 



— (T 
Ti 



0,H) 



K(T -> 0, H ) ~ J dr 2 (N^ orm ) /ttR 2 h - H 



(7) 



dr 2 (NJi orm ) 2 /TrR 2 H - H (6) 



where the contribution of M(tS) term is again neglected since 
its contribution is significant only near T c and becomes neg- 
ligible at low temperatures. The above field dependencies are 
for the ideally pure cases and we will show below, with nu- 
merical calculations, how these generic field dependencies are 
modified when the impurity scattering effect is included. 

Numerical results and discussions.- We calculate the 
l/Tx(T,H) and K(T, H) including the impurity scattering 
effect using the SCTA lflol, [Till . Once we calculate the im- 
purity induced selfenergy - normal and anomalous - correc- 
tions, £° mp (u;, r) and Sj mp (a;, r), respectively, we renormal- 
ize everywhere in the above formalism (Eqs. (1) -(3)) by 
p = " + Z° imp (u>,r) and A = A + E^w.r) (Z} mp van- 
ishes in the d-wave case due to the gap symmetry). In this 
paper, we considered only the strong coupling (unitary limit) 
non-magnetic impurity. The extension to the ±s-wave state 
for the calculations of N(lu, r) and M(lu, r) with two bands 
is straightforward and referred to Ref.B3[]. 

Figures 1 and 2 are the calculation results of the d-wave SC 
state. Fig. 1(a) shows the normalized l/Ti(T) of the pure 
d-wave state with various field strengths. First, the zero field 
data displays an exact T 3 power law for all temperatures be- 
low T c (for this we chose 2Ao/T c = 4). Applying fields, 
the low temperature behavior immediately changes to the T- 
linear power indicating that the zero energy DOS is created by 
the magnetic field. Increasing the field strength, the magni- 
tude of l/Ti(T) at low temperatures progressively increases. 
In the inset, we plotted l/Ti(T = 0.05T C ) vs JT/Jf c . As 
predicted in Eq.(4), the numerical data points (red circles) are 
perfectly fitted with H log H curve (black solid line; the fit- 
ting function is given in the inset with h = H/H C 2). 

Figure 1(b) shows the same calculations as in Fig. 1(a) but 
with impurity scattering. We chose the impurity concentra- 
tion parameter Y = ni mp /irN t ot(0) = 0.03Ao of the non- 
magnetic unitary scatterers. At zero field (black squares), 
1 /Ti (T) shows the typical T-linear at low temperatures due 
to the resonant (zero energy) impurity band but maintains the 
T 3 power law at higher temperatures. Applying fields, the 
overall line shape doesn't change but only the low temper- 
ature T-linear part progressively moves upward. It is clear 
that the zero energy DOS induced by the impurity scatter- 
ing is additive to the field induced zero energy DOS by the 
Volovik effect. This result was not a priori expected be- 
cause we calculated the position dependent impurity selfen- 
ergy S" mp (w, r) treating the Volovik effect and the impurity 
scattering on equal footing before we take the areal average. 
The plot of l/Ti(T = 0.05T C ) vs H/H c in the inset starts 
with a finite intercept at zero field, but still the data points are 
well fitted with aH log H + const, (black solid line). 

Figure 2 is the calculation results of the normalized Knight 
shift K(T) in parallel with the l/Ti(T) in Fig.l. It reflects 
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FIG. 1: (Color online) Normalized spin-lattice relaxation rates 
[l/Ti(T)/l/Ti(T c )] of the d-wave SC state with 2A /T C = 4 for 
various magnetic field strengths, H/H c2 = 0.0, 0.05, 0.1, 0.15, 0.2 
and 0.25, respectively, (a) The results without impurities (F = 0). 
(b) The results with the impurity concentration V — 0.03Ao of the 
unitary scatterer. Insets of each panel are the 1/21 vs H/H C 2 at a 
fixed low temperature T = 0.05T C . 
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FIG. 2: (Color online) Normalized Knight shift K(T)/K(T C ) of the 
d-wave SC state for various magnetic fields strengths, H/H C 2 = 
0.0, 0.05, 0.1, 0.15, 0.2 and 0.25, respectively, (a) The results with- 
out impurities (r = 0). (b) The results with the impurity concentra- 
tion r = 0.03Ao of the unitary scatterer. Insets of each panel are the 
K{H) /K N vs H/H c2 at a fixed low temperature T = 0.05T C . 



the same evolution of the DOS N(u, H) of the d-wave state 
with magnetic fields and impurities, as explained in Fig.l, 
onto K(T). Fig. 2(a) shows the normalized Knight shift K(T) 
of the pure d-wave gap with various field strengths. First, the 
zero field data displays the T-linear power law at low tem- 
peratures as expected for the nodal gap. Applying fields, the 
low temperature part immediately becomes flat because of 
the field induced zero energy DOS and progressively moves 
upward with increasing the field strength. The inset shows 
K(T = 0.05T C ) vs H/H c data and well fitted by the V# 
curve (black solid line) as shown in Eq.(5). Fig. 2(b) shows 
the same calculations of K(T) with impurity scattering as in 
Fig. 1(b). The overall behavior and evolution of it are now 
easily understood with the additive zero energy DOSs in- 
duced by impurities and by the magnetic field. The plot of 
K(T = 0.05T C ) vs H/H c data in the inset is well fitted by 
the ay/H + const, curve 11211 . 

Figures 3 and 4 are the calculation results of 1/Ti and K 
of the ±s-wave SC state. We assumed two gaps with differ- 
ent sizes, Ag/Ai = 0.4. Other parameters are: (1) A rather 
large value of 2Al/T c = 10 was chosen to produce the much 
steep drop of l/Ti(T) below T c ; (2) N L {0)/N s (0) = 2 is 
freely chosen although there is a correlation between two ra- 
tios, As/ Al and N L (0)/N s (0) As in the case of d-wave 
calculations, we considered only the unitary scatterers and as- 
sumed an equal strength between the interband and intraband 
scattering channels ll ill . 

Figure 3(a) shows the normalized l/Ti(T) of the pure ±s- 
wave state for various field strengths. First, the zero field 
l/Ti(T) (black squares) displays the exponential drop for 
all temperature range without the coherence peak. The ini- 
tial drop can be fitted by the T 6 power, much steeper -due 
to the choice of 2A L /T C = 10 - than the T 3 power that 
was often observed with many of Fe pnictide SC compounds. 
Applying fields, the low temperature behavior immediately 
changes to the quasi T-linear because of the zero energy DOS 
induced by the vortices as in the case of the d-wave gap. In- 
creasing the field strength, the magnitude of 1 jT\ (T) at low 
temperatures progressively increases. In the inset, we plotted 
l/Ti(T = 0.05T C ) vs H/H c and the numerical data points 
(red circles) are perfectly fitted with the iTlinear line (black 
dotted line) before it saturates when the small gap band be- 
comes completely collapsed beyond the critical field strength 
H*/H c2 ~ (A S /A L ) 2 as predicted in Ref.fi. 

Figure 3(b) shows the same calculations of 1/T\{T) but 
with impurity scattering. The value of the impurity con- 
centration r = 0.08Al was slightly larger than the critical 
concentration of the ±s-wave model UM (T* ~ 0.045Ai 
in this case), which would produce a perfect V-shape DOS 
and hence the T 3 power law at low temperatures. Therefore, 
in our case, the zero field data of 1/T\(T) (black squares) 
shows some power law in between T 3 and T-linear at low 
temperatures due to the combination of the finite zero energy 
DOS plus the V-shape DOS. Applying fields, the low temper- 
ature behavior immediately changes to the quasi T-linear as in 
the pure case of Fig. 3(a) and the overall magnitude progres- 
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FIG. 3: (Color online) Normalized spin-lattice relaxation rates 
[1 /Ti (T)/1/Ti (T c )] of the ±s-wave SC state for various magnetic 
fields H/H c2 = 0.0, 0.04, 0.08, 0.12, 0.16 and 0.2, respectively, (a) 
The results without impurities (T = 0). (b) The results with the 
impurity concentration T — 0.08Al of the unitary scatterer. Insets 
of each panel are the 1/Ti vs H/H C 2 at a fixed low temperature 
T = 0.05T C . Other parameters: 2A L /T C = 10, A S /A L = 0.4. 



FIG. 4: (Color online) Normalized Knight shift K(T)/K(T C ) 
of the ±s-wave SC state for various magnetic fields H/H C 2 = 
0.0, 0.04, 0.08, 0.12, 0.16 and 0.2, respectively, (a) The results with- 
out impurities (T = 0). (b) The results with the impurity concentra- 
tion r = 0.08Al of the unitary scatterer. Insets of each panel are 
the K{H) /K N vs H/H c2 at a fixed low temperature T = 0.05T C . 



sively increases with increasing the field strength. The plot of 
l/Ti(T = 0.05T C ) vs H/H c in the inset of Fig. 3(b) is much 
smoothened and there is no sharp kink at the saturation field 
H* fa Q.16H C 2\ but the i7-linear behavior still survives in the 
reduced region of low fields. 

Figure 4 is the calculation results of the normalized Knight 
shift K(T) of the ±s-wave SC state. Fig.4(a) shows the K(T) 
of the pure ±s-wave gap with various field strengths. First, the 
zero field data displays the exponentially flat behavior at low 
temperatures reflecting the full gap superconductor. Applying 
fields, the low temperature part moves upward and also de- 
velops a T-linear part up to Hj H C 2 ps 0.12. This behavior of 
K(T) ~ a+bT at low temperatures is the result of the combi- 
nation of the finite zero energy DOS plus the V-shape DOS - 
this is the typical feature of the multi-band ±s-wave SC state 
with the resonant impurity scattering! 1 1] but it is now gener- 
ated by the magnetic fields. The inset shows K(T = 0.05T C ) 
vs H/ H c data which is well fitted by the i7-linear line (black 
dotted line) as shown in Eq.(7) at low fields. 

Figure 4(b) shows the same calculations of K(T) includ- 
ing impurity scattering of T = 0.08Al. The overall behavior 
and evolution of it are again easily understood with the ad- 
dition of the impurity induced DOS and the magnetic field 
induced DOS. The inset plot of K(T = 0.05T C ) vs H/H c 
data is much rounded and no kink feature is visible as in the 
inset of Fig. 3(b). In general, all data of K(T) calculations 
in Fig. 2 and Fig. 4 faithfully reflect the evolution of the total 
DOS Ntot(w) at low frequencies by magnetic fields and im- 



purity scattering as can be seen from Eq.(3). 

Conclusion. - In summary, we extended the study of the 
Volovik effect on the NMR properties both in the d-wave and 
±s-wave gap states. We derived the generic power laws of the 
field dependencies of -> 0,H) and K(T -> 0,H) 

in pure cases. We also numerically calculated the full temper- 
ature dependencies of l/Xi(T, H) and K (T, H) for various 
field strengths, for both SC states, including impurity scatter- 
ing. Main findings are: (1) the field dependence power laws 
of 1/T X (T -> 0,H) and K(T -> 0,H) are practically un- 
changed with impurity scattering in the d-wave case except 
for a constant shifting ll 1211 . In the ±s-wave case, while the 
changes of the overall field dependencies are more significant, 
the generic power laws survive in a good part of low field re- 
gion; (2) the Volovik effect acts as an equivalent pair breaker 
as the unitary impurity scatterers to creates the zero energy 
DOS and hence to modify the low temperature power laws of 
l/Ti(T) and K(T). It implies that the Volovik effect always 
needs to be taken into account for the analysis of the NMR 
experiments. 
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